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The purpose of this mock data and science challenge is to prepare the data analysis and science 
interpretation for the second generation of gravitational-wave experiments Advanced LIGO-Virgo 
in the search for a stochastic gravitational-wave background signal of astrophysical origin. Here 
we present a series of signal and data challenges, with increasing complexity, whose aim is to test 
the ability of current data analysis pipelines at detecting an astrophysically produced gravitational- 
wave background, test parameter estimation methods and interpret the results. We introduce the 
production of these mock data sets that includes a realistic observing scenario data set where we 
account for different sensitivities of the advanced detectors as they are continuously upgraded toward 
their design sensitivity. After analysing these with the standard isotropic cross-correlation pipeline 
we find that we are able to recover the injected gravitational-wave background energy density to 
within 2(7 for all of the data sets and present the results from the parameter estimation. The results 
from this mock data and science challenge show that advanced LIGO and Virgo will be ready and 
able to make a detection of an astrophysical gravitational-wave background within a few years of 
operations of the advanced detectors, given a high enough rate of compact binary coalescing events. 


I. INTRODUCTION 

According to various cosmological and astrophysical 
scenarios, we are bathed in a stochastic gravitational 
wave background (SGWB). Proposed theoretical cosmo¬ 
logical models include the amplification of vacuum fluc¬ 
tuations during inflation[TH3], pre Big Bang models [lHS]j 
cosmic (super)strings [71-fTU]. and phase transitions nu¬ 
lls]. In addition to the cosmological background [nnm, 
an astrophysical contribution m is expected to result 
from the superposition of a large number of unresolved 
sources, such as core collapse supernovae to neutron stars 
or black holes [I2H2 q], rotating neutron stars [2TJ[22|, in¬ 
cluding magnetars , phase transitions m , or ini¬ 

tial instabilities in young neutron stars PHH5D] . or com¬ 
pact binary mergers pTIVUIfi) . 

Many of these models are within reach of the next gen¬ 
eration of gravitational-wave (GW) detectors such as Ad¬ 
vanced LIGO Hanford (H) and Livingston (L) [36| , which 
are expected to start collecting data in 2015, and Ad¬ 
vanced Virgo (V) [37|, which will begin collecting data 
in 2016. These detectors are expected to have a final 
sensitivity 10 times better than that of the initial de¬ 
tectors [351 [3S| 5 which will be achieved over a period of 
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several years of continued upgrades. The detection of a 
cosmological background would provide very important 
constraints on the first instant of the Universe, up to 
the limits of the Planck era and the Big Bang, while the 
detection of an astrophysical background would provide 
crucial information on the physical properties of the re¬ 
spective astrophysical populations, the evolution of these 
objects with redshift, the star formation history or the 
metallicity [THl . 

In order to prepare and test our ability at detecting 
the SGWB and interpreting valuable information from 
the data, we are conducting a series of mock data and 
science challenges (MDSC), with increasing degrees of 
complexity. In this first paper, we focus on a SGWB cre¬ 
ated by all the unresolvable compact binary coalescences 
(GBC) such as binary neutron stars (BNS), neutron star- 
black holes (NSBH) or binary black holes (BBH), up to 
a redshift z = 10, which may dominate within our search 
frequency range. Such a background may have a realistic 
chance of being detected after a few years of operation 
of the advanced detectors (33). The observation of the 
SGWB will complement individual detections of a few to 
a few tens of CBG events per year [33] at close distances 
up to a few hundred Mpc. 

For this study we produced multiple year-long data 
sets in the form of time series for the three advanced 
LIGO/Virgo detectors, containing both instrumental 
noise and the GW signal from a large number (10^ — 10^ 
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per year) of compact binary sources out to redshift 
z = 10. This is done using a data generation package 
that was initially developed for the Einstein Telescope 
MDSC [inillT]. The data sets are then analysed using a 
cross-correlation (CC) analysis pipeline in order to mea¬ 
sure the total GW energy density, flgw)/), of all the GW 
signals that make up the stochastic background [?5] . We 
then use the results from these analyses to perform pa¬ 
rameter estimation to try to determine some of the pa¬ 
rameters of the injected populations, such as the average 
mass of all the sources, or the rate at which these binaries 
coalesce. 

This paper is organised into the following sections. In 
SectionHH we introduce the mock data sets and the meth¬ 
ods by which we generate them. In SectionjTTlj we briefly 
describe the analysis methods used to detect the stochas¬ 
tic signal. In Section |IVj we discuss how we use the re¬ 
sults from the analysis to estimate various astrophysical 
parameters. In Section [Vj we present the results from 
various analysis runs of the different mock data sets. Fi¬ 
nally in Section [VI[ we present our conclusions. 

II. MOCK DATA 

In this section, we introduce the mock data sets that 
we will be analysing as part of this investigation, as well 
as the data generation program we use to produce them. 
Initially we explain the various steps that are used to 
produce the mock data before detailing each of the data 
sets that will be considered as part of this MDSC. Finally 
we show how one can consider an astrophysical SGWB 
as the superposition of many unresolvable single sources. 
The main data sets are split into two subsets, the first 
being produced with Gaussian noise, and the second pro¬ 
duced using “glitchy” data taken from the initial LIGO 
S5 and initial Virgo VSRl science runs re-coloured to 
have the sensitivity of the advanced detectors [JS]. The 
GW signals injected into each set are the same for both 
the Gaussian and re-coloured noise, allowing us to make 
a direct comparison of how the analysis pipeline will be¬ 
have in an ideal and in a more realistic case. 


A. Mock Data Generation 

The mock data generation package used here was orig¬ 
inally developed for the Einstein Telescope MDSC [46] 
where one would expect to be able to make detections 
of individual sources out to z « 3.8 for BNS and even 
further for higher mass systems such as NSBH or BBH. 
Being able to realistically represent the population of 
sources at high redshift is essential when considering an 
astrophysical SGWB signal for the advanced detectors. 
This is because we expect very few CBC events to be 
directly detectable but the large number of unresolvable 
sources, when considering the whole universe, will all con¬ 
tribute to the SGWB. This is clearly seen in Fig. [^where 



FIG. 1: Source redshift probability distribution for BNS 
(blue) and BBH (red) where different delay times between 
the formation and merger of the binary system are consid¬ 
ered. We also plot the horizon distance to BNS (dashed blue) 
and BBH (dashed red) for the advanced detectors with their 
design sensitivity which are defined as 445 Mpc (z ~ 0.1) and 
2187 Mpc (z ~ 0.4) respectively [45] . 

we show the redhsift probability distribution, which is ex¬ 
plained in the next section, of BNS (blue) and BBH (red) 
which take into account different delay times between the 
formation and merger of the inspiralling systems. We also 
plot the maximum horizon distance of the advanced de¬ 
tectors to directly detect individual BNS (blue dashed) 
and BBH (red dashed) signals [50] . We now describe how 
we generate and add the large number of GW signals to 
the detector data streams using both Gaussian and re¬ 
coloured noise. 


1. Simulation of GW CBC signals 

The Monte Carlo procedure we use to generate a pop¬ 
ulation of compact binaries up to redshift of z = 10 is 
described in detail in [461152] . Here we summarize the 
main steps of the simulations. 

The coalescence of a compact binary occurs after two 
massive stars in a binary system have collapsed to form a 
neutron star or a black hole^ and have inspiralled through 
the emission of gravitational waves. 

The coalescence rate in Mpc“^ yr“^ is given by [53U56] 

Pciz,td) = X [ Pitd)dtd. (1) 

J 1 + z/ 

In this expression, is the star formation rate (SFR), 
measured in Mq Mpc“^ and the factor l/{l-\-Zf) ac- 


^ We neglect the possible production of compact binaries through 
interactions in dense star systems. 
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FIG. 2: Left — Evolution of the Advanced LIGO sensitivities over the early (blue), middle (red), late (green) and design (black) 
phases m- Right — Evolution of the Advanced Virgo sensitivities over the early (blue), middle (red), late (green) and design 
(black) phases |51 |. 


counts for time dilation due to the cosmic expansion. The 
redshift 2 : describes when our compact binary merges, 
2 / is the redshift at which its progenitor massive binary 
formed, and P{td) is the probability distribution of the 
delay between 2 and Zf, which is the sum of the time 
from initial binary formation to evolution into a com¬ 
pact binary, plus the merger time by emission of gravi¬ 
tational waves. The parameter A (in Mq^) is the mass 
fraction that goes into the formation of the progenitors 
of compact binaries. The local coalescence rate at 2 = 0, 
/Oc(0,td)) is one of the parameters of our simulations and 
is denoted by po- 

The merger rate in the redshift interval [ 2 , 2 - 1 - d 2 ] is 
obtained by multiplying pdz, td) with the element of co¬ 
moving volume 


Hubble parameter Hq — 70 km s ^ Mpc Follow¬ 
ing [31I1T] we assume a distribution of the delay of the 
form Pitd) oc 1/td with a minimal delay of 20 Myr for 
BNS and 100 Myr for BBH, as suggested by population 
synthesis [55H5D] . 

We then proceed as follows for each source: 

• The arrival time of each GW signal is selected 
from a Poisson distribution. Here, the difference in 
arrival time, t = where k is the current 

event, is drawn from the exponential distribution 
P(t) = exp(— t/t) where r is the average time be¬ 
tween successive events. The average waiting time 
between signals is calculated by taking the inverse 
of the coalescence rate, Eq. ([^, integrated over all 
redshifts 


where 


di? , N . / N dH , . 
— (2,td) = pc{z,td)-^{z}, 


dV c r{zd 

-^( 2 ) -47r—^(^^ 2 )’ 


( 2 ) 


(3) 


where c is the speed of light in a vacuum, Hq is the 
Hubble constant and 


Jo 


dz' 
0 


(4) 


where 


E[fl,z) — •\/Ha + Hm,(l + 2 )^ . (5) 

Here, we select the SFR given in m and use a stan¬ 
dard ACDM cosmology with = 0.3, Ha = 0.7 and 



di? 

d^ 


iz,td) dz 



( 6 ) 


• The redshift at the point of coalescence, 2 , is se¬ 
lected from a probability distribution p{z, td) con¬ 
structed by normalizing the coalescence rate in the 
interval [0,10] (see Fig. 


p{z,td) = T — {z,td). (7) 

• The SGWB analysis is not sensitive to the width of 
the distribution of the masses, only to the average 
chirp mass, JA, of the system. This is a combi¬ 
nation of the two component masses, mi and m 2 , 
given by 


(mim2)3/5 
(M)i/5 ’ 


( 8 ) 
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where M = mi +m 2 is the total mass of the system. 
Because of this, we choose a single value for the 
component masses for each of the systems being 
considered: 1.4M0 for neutron stars and IOMq for 
black holes. 

• The sky position, Cl, is selected from an isotropic 
distribution across the whole sky. The cosine of 
the inclination angle of the orbital plane to our line 
of sight, i, the GW polarisation angle, ip, and the 
phase angle at the time of coalescence, (po, are all 
chosen from uniform distributions. 


and hx are given in [531. The total SNR for the 
GW detector network is then given by 

P^='Epa, (13) 

A 

where A is the sum over all detectors in the net¬ 
work. Any signals that pass a network threshold 
SNR pt are then ignored by the SGWB search, 
where we set the network SNR threshold to 12. 


• We next calculate the theoretical signal-to-noise ra¬ 
tio (SNR), p, of the inspiral signal to determine if 
it is individually detectable by the standard LIGO- 
Virgo CBG search pipeline jUj [551 |5T| (52]. The 
SNR produced by the inspiral phase of the wave¬ 
form for GBCs is given by 


p2=4 


F+{Cl,iP)h+{f)+Fx{Cl,iP)hxif) 

Snjf) 


d/, 

(9) 


where A+ and Fx are the antenna response func¬ 
tions to the two GW polarisations originating from 
sky position and with polarisation angle ip, and 
and hx are the signal amplitudes in the fre¬ 
quency domain for the two polarisations. In the 
Newtonian regime, before the last stable circular 
orbit, h+ and hx are given by 


• Finally, for any surviving sub-threshold events, we 
produce the waveforms that are then added to the 
detector data steams. Here, we have chosen to 
use the TaylorT4 time-domain waveform up to 3.5 
post-Newtonian order in phase, and the most dom¬ 
inant post-Newtonian lowest order for amplitude, 
for the BNS signals. For the case of BBH signals, 
we choose the EOBNRv2 waveform produced from 
numerical relativity, which includes the merger and 
ringdown of the two coalescing black holes. This is 
up to 4th post-Newtonian order for phase and the 
lowest order for amplitude |5i] . 

Once the time series data has been produced containing 
all the sub-threshold injections, we add either Gaussian 
noise or re-coloured noise data to produce the final mock 
data sets. 


2. Simulation of Gaussian noise 


(1 + cos^ i) 7/6 

z ^ J , 


h+{f) 

hxif) = /iz cost/"'^/®. 


( 10 ) 

( 11 ) 


where 


hr = \ — 


24 TT'^/^c^Gdi^{z) 


( 12 ) 


In the above equations G is the gravitational con¬ 
stant, cIl is the luminosity distance to the source 
at redshift z, /min is the starting frequency, which 
we select to be 10 Hz, Sn is the detector’s noise 
power spectral density (PSD) (see Fig.[^, and /isco 
is the frequency of the last stable circular orbit, 

^3 


/isco — 


For BNS signals it is enough 


63/27rGM’ 

to consider the waveform up until this point as the 
SNR contribution of the inspiral phase is dominant. 
However, for BBH signals we must also consider the 
contribution from the merger and ringdown of the 
waveform. The modifications to the calculations of 


Because each of the detectors being considered in this 
MDSG are well-separated in space, we assume that there 
will be no correlated noise between any of them; so the 
noise is simulated independently for each of the detec¬ 
tors [551155] . We do this by generating a mean zero, unit 
variance Gaussian time series which is then Fourier trans¬ 
formed into the frequency domain. This is then coloured 
using the PSD of the detector sensitivity of either aLIGO 
or AdVirgo (see Fig. and is then finally Fourier trans¬ 
formed back into the time domain. To prevent any poten¬ 
tial discontinuities of the data between adjacent segments 
of data, we taper the noise curve away to zero at frequen¬ 
cies below lOHz and above 512Hz. This is then added to 
the time series containing the injected GW signals. 


3. Re-coloured noise 

To more accurately imitate the noise likely to be 
present in the advanced detectors, we also re-colour ini¬ 
tial S5 LIGO and Virgo VSRl data to have the sensi¬ 
tivity of aLIGO and AdVirgo. This has the benefit of 
preserving non-stationary noise artifacts while exploring 
the sensitivity of the pipelines to the signals. The noise 
spectra are estimated from the year of data and then 
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averaged. The data is then re-coloured with a transfer 
function corresponding to the advanced detector power 
spectra divided by this averaged spectrum. The same 
GW signals are then added to these re-coloured time- 
series. 


B. Mock Data Sets 

To ensure that we have a detectable signal and to re¬ 
duce computational costs, we have selected astrophysical 
models corresponding to the most optimistic rates given 
in [JS] instead of using longer observational times to ob¬ 
tain the same SNR. 

In total, we generate 5 data sets, of duration one year, 
all of which are produced with both Gaussian noise and 
re-coloured noise. For these data sets, we use the nomi¬ 
nal design sensitivities of Advanced LIGO and Virgo, but 
in order to account for the improvement of the sensitiv¬ 
ity, we also produce a sixth, observing scenario data set 
consisting of 3, 6, and 9 months and 1 and 3 years cor¬ 
responding to the early, middle, late and design phases 
of the advanced detectors (see Fig. [^. This observing 
scenario, consisting of 5.5 years of data, should cover the 
full advanced detector observing period from mid 2015 
until the end of 2022. Details of all the data sets are 
found in Table ||] and details on the rates and expected 
number of events are found in Table IE 

One part of the investigation is to see how having sta¬ 
tistically different sources effects our analysis. When the 
rate and the duration of the GW events are large, the 
sources overlap each other creating a GW signal contin¬ 
uous in time (there is always a source present) which 
is Gaussian in nature due to the central limit theorem. 
However, for smaller rates or shorter waveforms, the time 
interval between successive events increases resulting in 
a non-continuous and non-Gaussian signal [niinzi- Ana¬ 
lytical calculations made in [68] and for the Einstein Tele¬ 
scope MDSCs [461 m] suggest that when making a mea¬ 
surement of an astrophysical gravitational-wave back¬ 
ground one just needs to consider the total number of 
coalescing events, with their relevant signal amplitudes, 
that occur within the observational period. The nature of 
the signals themselves will have no effect on the estima¬ 
tion of the SGWB. This is explained in greater detail later 
in Section III Cl However these results have not been in¬ 
dependently verified with the use of simulated data. We 
now describe each of these data sets being considered in 
this investigation. 

1. Main data sets 

Data set 0 is our control test for both the Gaussian 
and re-coloured data sets. The data streams for each of 
the detectors contain no coincident signals so there will 
be no correlated signals between any of the detectors. 
Thus, the results from the analysis of data set 0 should 


give us an accurate measurement of the expected error 
bars for each of the following data sets. 

In data set 1, we have generated a large number of 
BNS signals with a merger rate of 10 coalescences per 
Mpc^ per Myr that are injected into both the Gaussian 
and the re-coloured noise. Any individual events that 
surpasses a network SNR threshold value, as described 
in Section |H A 1[ are removed as it is possible for these 
signals to bias the results of the analysis. But given the 
expected number and length of the waveform when com¬ 
pared to the overall length of the observing time this 
effect is negligible, as demonstrated later in this paper. 
The top plot of Fig. [^ shows a 1000s segment of the time 
series produced from BNS contributing to data set 1. It 
is clearly seen that the GW signal is continuous at all 
times, so with this data set we investigate the affect that 
a continuous SGWB signal will have on our analysis. 

Data set 2 contains the exact same sources with both 
Gaussian and re-coloured noise that are included in set 
1 except that here we have not removed the loudest in¬ 
dividually detectable events from the time series. This 
is to test by how much our results can be affected if we 
include the loud events in the analysis by comparing the 
results against that of the first data set. 

In data set 3, we generate a number of BBH signals 
using a merger rate of 0.3 coalescences per Mpc^ per 
Myr that are injected into both Gaussian and re-coloured 
data. With this data set we investigate the possible ef¬ 
fects that a non-continuous (popcorn) SGWB will have 
on our analysis. The middle plot of Fig. [^ shows a 1000s 
segment of the time series produced from BBH contribut¬ 
ing to data set 3. We see that, due to the shorter wave¬ 
form lengths and lower coalescence rate, the GW signals 
are non-continuous or more “popcorn” like. 

In data set 4 we have generated a number of both BNS 
and BBH signals, using the same rates stated above, 
which are injected into both Gaussian and re-coloured 
data. This set is to test the behaviour of the analysis and 
parameter estimate when analysing data from a SGWB 
having a contribution from more than one source. The 
bottom plot of Fig. [^ shows a 1000s segment of the time 
series produced from BNS and BBH signals contributing 
to data set 4. It is clearly seen that the GW signal is 
continuous at all times whilst still having louder popcorn 
like bursts from the BBH signals. 


2. Observing scenario 

The observing scenario data set is designed to realis¬ 
tically simulate the data that we would expect to obtain 
from the network of advanced detectors over the initial 
several years of operations. For this set, we generate a 
large number of BNS signals using a merger rate of 2 co¬ 
alescence per Mpc^ per Myr that are injected into Gaus¬ 
sian data. This is a lower rate than is used in set 1 as we 
are using a longer observational period. This Gaussian 
data differs from the previous data sets as we change the 
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PSDs for each of the detectors at different stages to rep¬ 
resent the improvements in sensitivities that are expected 
to be obtained in each observing run. Examples of this 
are shown in the left-hand plot for aLIGO, and right- 
hand plot for AdVirgo of Fig. In reality we should 
take the final sensitivity of one phase as the initial sen¬ 
sitivity of the next, which would then gradually decrease 
to the next final sensitivity. However, here we consider 
the ideal case of taking the final sensitivity for the full 
duration of each phase. These observing runs are set out 
as follows [5T] : 

1. 01, 2015: This will consist of a 3-month obser¬ 
vational run with both LIGO detectors (HL) with 
Early aLIGO sensitivity. 

2. 02, 2016-17: This will consist of a 6-month obser¬ 
vational run with all three detectors (HLV), where 
the HL detectors will have Mid aLIGO sensitivity 
and V will have Early AdVirgo sensitivity. 

3. 03, 2017-18: This will consist of a 9-month obser¬ 
vational run with all three detectors (HLV), where 
the HL detectors will have Late aLIGO sensitivity 
and V will have Mid AdVirgo sensitivity. 

4. 04, 2019: This will consist of a year long obser¬ 
vational run with all three detectors (HLV), where 
the HL detectors will have the final design aLIGO 
sensitivity and V will have Late AdVirgo sensitiv¬ 
ity. 

5. 04 (continued), 2020-22: This will consist of a 3 
year long observational run with all three detectors 
(HLV), where the HL detectors will have the final 
design aLIGO sensitivity and V will have the final 
design AdVirgo sensitivity. 


C. Astrophysical stochastic background from CBC 

The spectrum of the SGWB is usually characterized 
by the dimensionless parameter 


f^gw(/) 


1 dpgw 
Pc d In / ’ 


(14) 


where pgw is the gravitational-wave energy density, / the 

3^2 

frequency in the observer’s frame and, pc = , is the 

critical energy density needed to close the Universe today. 
The GW spectrum from the population of extra-galactic 
compact binaries is given by the expression 


f^gw(/) = —fFif), 

PcC 


(15) 
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FIG. 3: Top — A 1000s segment of the time series for BNS 
signals using the higher LIGO rate of 10 Mpc“® Myr“^. This 
signal appears continuous. Middle — A 1000s segment of the 
time series for BBH signals of mass 10 -I- 10 Mq, using the 
higher LIGO rate of 0.3 Mpc“^ Myr“^. This signal appears 
non-continuous (popcorn). Bottom — A 1000s segment of the 
time series for mixed BNS and BBH signals with higher LIGO 
rates and BH mass of IOMq . This signal appears as a popcorn 
background from BBH on top of a continuous background 
from BNS. 


where F(f) is the total flux. The total flux is the sum of 
the individual contributions 
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TABLE I: Table describing the data sets that are produced as part of the MDSC. The first column is the reference number 
of the data set. The second column indicates if the data set is produced with just Gaussian noise or with re-coloured noise as 
well. The third column shows what sources are injected into the data set. The forth column gives the rate of events used, see 
Table [III The fifth column gives the length of the data set. 


Data Set 

Noise 

Sources 

Rate 

Tobs 

0 

Gaussian & Re-coloured 

- 

- 

1 year 

1 

Gaussian & Re-coloured 

BNS (sub-threshold) 

10 Mpc“^Myr“^ 

1 year 

2 

Gaussian & Re-coloured 

BNS (all) 

10 Mpc“®Myr“^ 

1 year 

3 

Gaussian & Re-coloured 

BBH (sub-threshold) 

0.3 Mpc“^Myr“^ 

1 year 

4 

Gaussian & Re-coloured BNS -h BBH (sub-threshold) 10+0.3 Mpc ^Myr ^ 

1 year 

Observing Scenario 

Gaussian 

BNS (sub-threshold) 

2 Mpc ^Myr ^ 

5.5 years 


TABLE 11: A list of compact binary coalescence rate den¬ 
sities as given in [IS]. The first column labels whether a 
merger rate is optimistic (i?high), realistic (i7reaiistic) or pes¬ 
simistic (i?iow). The second column gives the rates of coalesc¬ 
ing events. The third column gives the average time between 
successive events. The hnal column gives the total number of 
events out to z = 10 that are expected to occur per year. 


(16) 


Expected Rate po (Mpc ^ Myr ^) 

r (s) 

Nevents (yf ) 

BNS 

-Rhigh 

10 

1.35 

2.3 X 10^ 

-Rrcalistic 

1 

13.5 

2.3 X 10® 

-Rmedium-low 

0.1 

135 

2.3 X 10® 

Rlow 

0.01 

1350 

2.3 X 10"^ 

BBH 

-Rhigh 

0.3 

64.7 

4.9 X 10® 

-Rrcalistic 

0.005 

3880 

8133 

-Rmedium-low 

0.001 

19400 

1627 

-Rlow 

0.0001 

194000 

163 


where N is the total number of coalescences in the data 
and k is the index of the individual coalescence. The 
normalization factor assures that the flux has the 
correct dimension, Tobs = 1 yr being the length of the 
data sample. 


In the Newtonian regime, before the last stable circu¬ 
lar orbit, the Fourier transforms and hx are given by 
Eqs. (10 - 12). This gives for the energy density param¬ 
eter [68j 


^^gw(/) 


5^2/3g5/3^5/3 ^ + 


(1 -I- COS^ Lk)"^ 


-I- COS^ 



(17) 


This equation is valid for BNS signals where we have 
considered only the inspiral phase, but for BBH signals, 
there is an extra contribution coming from the merger 
and ringdown phases. Fig. 13 shows f2gw(/) for the pop¬ 
ulation of BNS (blue), and BBH (red) in the mock data 
sets described in the previous section. The plot in black 
corresponds to the sum of the signal from BNS and BBH. 
For BNS and BBH signals, Hgw(/) increases as from 
the inspiral phase (then as from the merger phase for 
BBH) before it reaches a maximum and decreases dra¬ 
matically. The peaks occur at frequencies corresponding 
roughly to the f\sco and the end of the ring-down phase 
at 0 ^ 1.5 where the coalescence rate is maximal. The 
amplitude of the background scales with both the rate of 


coalescing events and the average chirp mass of all the 
signals. It is larger for the BNS background (data set 
1) than for the BBH contribution (data set 3) because, 
even though the chirp mass is smaller, the rate we have 
considered is larger. 


III. SGWB SEARCH 

In this section, we briefly describe the cross-correlation 
(CC) method by which we analyse the data (please see 
references such as [40l|48l|69] for a complete treatment). 
The optimal strategy to search for a Gaussian (or contin¬ 
uous) SGWB is to cross-correlate measurements of mul- 

























FIG. 4: Left — Energy density Qg^{f) for the population of BNS (blue), BBH (red), and the combination of the two populations 
(black) as used in data sets 1, 3 and 4 respectively. The plots are calculated using the higher LIGO rates of 10 Mpc“® Myr“^ 
(BNS) and 0.3 Mpc“^ Myr“^ (BBH). We also plot the reference frequency for which we report all the results in this investigation, 
shown by the black dashed line. Right — Fraction of the total theoretical SNR for the three d etec tor pairs, HL (blu), HV (red), 
and LV (green). This is calculated using the frequency band [10 — 250]Hz, as given by Eq. (22l, and assuming which is 
used for the analysis. We also plot the reference frequency for which we report all the results in this investigation, shown by 
the black dashed line. 


tiple detectors, s'i(/) and S 2 (/). When the background 
is assumed to be isotropic, unpolarized and stationary, 
the cross correlation product is given by [53] 


/ OO 

Sl*(l/l)S2(|/|)Q(/)d/, (18) 

-OO 

and the expectation value of Y is 


SNR = 



3^0^ 

I V 


IOtt^ 


obs 




ppiiumm 


df 


( 22 ) 


where we usually assume a filter of the form 


^OO -I 

< ^ ^ ^f^gw(l/l)7(l/l)Q(/)d/, 

(19) 

where 


Qif) 


\f\^Pi{\f\)P2{\f\y 


( 20 ) 


is the optimal filter that maximizes the SNR, is 

the energy density in GW as defined in Eq. ([M]), si and 


s ~2 are the detector output streams from both detectors in 
the frequency domain, Pi{f) and P 2 {f) are the detector 
power spectral densities of the two detectors and 7 (/) is 
the normalized overlap reduction function m, charac¬ 
terizing the loss of sensitivity due to the separation and 
the relative orientation of the detectors. 

The expected variance, which is dominated by the 
noise, is given by 


rp POO 

Pi(|/|)P2(|/|)|Q(/)pd/, (21) 

and the expected (power) SNR of the CC statistic in the 
frequency range /min — /maxi for an integration time Tobs 
is given by m 


f^gw(/) = ^M/freiT- (23) 

For this MDSC, we set /min = lOHz, /max = 250Hz, 
use a reference frequency of /ref = lOOHz and set 0 = 2/3 
as this is the theoretical value produced from the inspiral 
phase of CBCs. In the right hand plot of Fig. [^ we show 
the fractional SNR build-up as a function of frequency for 
different detector pairs. We see that for all three pairs we 
reach to nearly 100% of the total SNR by 120Hz, which 
is contained well within the limits we have set above for 
the analysis. 


IV. PARAMETER ESTIMATION 


Parameter estimation of signal models requires, at 
first, GW detection with high significance. In this analy¬ 
sis, we use a method for parameter estimation of a SGWB 
background. We seek to address the question of how 
well we can fit the model parameters. As a concrete 
example, we show the recovery of parameters from the 
MDSC injection sets using both a power-law, a, and CBC 
model. We show how to estimate parameters such as a 
SGWB amplitude and the CBC coalescence rate. To do 
so, we use a method presented in [33], that introduced 
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a maximum likelihood technique to simultaneously es¬ 
timate multiple parameters of SGWB models using CC 
data from pairs of GW detectors. This technique was 
used on recent results from LIGO to produce the first 
simultaneous limits on multiple parameters for power- 
law and CBG models of the SGWB, and to estimate the 
sensitivity of second-generation GW detectors to these 
models. 

The likelihood function is defined as 


L{Yi,a\9) oc 


1 

■9 ^ 


21 


(24) 


where \ 0) is the template spectrum that we are try¬ 

ing to fit by varying the parameters 0, the sum runs over 
frequency bins /i, which we set to be 0.25Hz, and Yi and 
af are the estimator and variance in the frequency bin i. 
The two methods, one a stochastic template-based anal¬ 
ysis and the other a GG analysis, are very similar. Tradi¬ 
tional SGWB cross-correlation searches have adopted a 
particular power-law model, assuming a specific spectral 
index value and searching over the spectral amplitude. 
The template-based maximum-likelihood estimator in¬ 
stead generically incorporates any model for the purposes 
of both detection and parameter estimation. Therefore, 
it may be particularly useful for compact binary back¬ 
grounds. 


mass fraction parameter, proportional to the local CBG 
rate per unit volume. 


V. RESULTS 

We now present the results from our analyses of the 
various mock data sets. We first discuss the results from 
the analysis of the Gaussian noise data (data set 0), 
which we consider to be an ideal case. We compare these 
against the results from the re-coloured noise which can 
be considered as a more realistic case. We then detail 
the results from the simulation of the observing scenario. 
Finally we show the results from the parameter estima¬ 
tion. 

The results of the Gaussian and re-coloured data sets, 
where we use the same set of injections for both, are 
reported in Table m The first column lists the three 
detector pairs as well as the combined results which is 
the weighted sum of the three pairs where the combined 
point estimate is calculated using 


^combined — 


Yabc^ab 

AB 


AB 


r-2 

^AB 


(25) 


where AB run over the three possible detector pairs and 
the combined error is given by 


A. SGWB Models 

The amplitude and the frequency dependence of the 
SGWB spectrum depend on the physics of the model 
that generated it. For example, in the CBG model, the 
spectrum is determined by the coalescence rate of binary 
systems throughout the universe and by the distribution 
of their chirp masses. Past SGWB searches, performed 
using data from the initial LIGO and Virgo detectors [101 
1151175] , assumed a power-law model (see Eq. (231) and set 
limits only on the amplitude flref for fixed values of the 
spectral index a and of the reference frequency /ref. This 
is reasonable as most SGWB models predict a power-law 
dependence in the LIGO-Virgo frequency band. 

As discussed above, compact binary coalescences are 
among the most promising sources of gravitational waves 
for ground-based gravitational-wave detectors. While de¬ 
tections of individual compact binaries are possible, an¬ 
other possibility is the detection of contributions from all 
GBCs in the universe to a SGWB. The model we use is 
the average version of Eq. ( |17| ) where the discrete sum 
over the sources is replaced by an integral over the red- 
shift, masses, sky position, inclination angle and ip. It 
was shown in [44] that it was sufficient to use only the 
average chirp mass Ai to determine the spectrum. The 
model then consists of A4 and A, which we recall is the 


-2 _ 

^combined / . 

AB 


AB- 


(26) 


The second, third and fourth columns give the estimated 
GW energy density, the error on this estimate and the 
SNR of the measurement for the data sets using Gaus¬ 
sian noise. The fifth, six and seventh columns give the 
corresponding results for the data sets using re-coloured 
noise. We define an SGWB signal as being detectable 
once it passes an SNR threshold value of 3, as defined in 
Eq. 


A. Gaussian 


The results from the analysis of all the Gaussian data 


sets are presented in the left-hand section of Table III 


The first results we highlight are from data set 0, which 
is our control set as it consists of just independent Gaus¬ 
sian noise. The measurement of flo, for the three detector 
pairs and the combined result give very low estimates and 
are the results of statistical fluctuations in the Gaussian 
data. The main result from this data set is the measure¬ 
ment of the error which we note is consistent with the 
error measurements of each of the other Gaussian noise 
data sets. 

Data sets 1 and 2, where we only consider a popula¬ 
tion of BNS with a high merger rate, were designed to 
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TABLE III: Results from all mock data sets. The first column indicates the detector pair used in the analysis. The second 
column gives the estimated value of from the Gaussian data sets. The third column gives the error on the measurement 
from the Gaussian data sets. The fourth column gives the SNR of the detection from the Gaussian data sets. The fifth column 
gives the estimated value of Qa from the re-coloured data sets. The sixth column gives the error on the measurement from 
the re-coloured data sets. These are all nearly identical across the data sets because the contributions of signal to the overall 
noise background are minimal (of order 0.1% when compared to the instrumental noise). This was verified by performing the 
same analysis on a small subset of the data but containing only signal. The final column gives the SNR of the detection from 
the re-coloured data sets. The table is divided by horizontal rows for the various data set with the injected value of fla, as 
calculated by Eq. (171, also given. 


Detector Pair 

Point Estimate Error 

SNR 

Point Estimate Error SNR 


Gaussian 

Re-coloured Noise 


Set 0 (Noise only): Qa = 0 


HL 

-8.093 X 10“"“ 

1.473 X 10““ 

-0.55 

-1.119 X 10““ 

1.683 X 10““ 

-0.66 

HV 

-1.04 X 10“® 

1.139 X 10“® 

-0.91 

-1.12 X 10“® 

1.407 X 10““ 

-0.8 

LV 

-1.47 X 10“® 

1.042 X 10“® 

-0.14 

-1.765 X 10““ 

1.38 X 10“® 

-0.13 

Combined 

-9.769 X 10“"“ 

1.447 X 10““ 

-0.68 

-1.268 X 10““ 

1.659 X 10““ 

-0.76 


Set 1 (BNS): 

= 1.364 X 10“® 



HL 

1.512 X 10“® 

1.474 X 10“® 

10.26 

1.455 X 10““ 

1.683 X 10““ 

8.65 

HV 

7.706 X 10“'^ 

1.139 X 10“® 

0.68 

-9.858 X 10““ 

1.235 X 10““ 

-0.8 

LV 

5.491 X 10“® 

1.042 X 10“® 

0.53 

7.451 X 10““ 

1.235 X 10““ 

0.6 

Combined 

1.481 X 10““ 

1.448 X 10““ 

10.23 

1.399 X 10““ 

1.653 X 10““ 

8.46 


Set 2 (BNS): 

= 1.411 X 10““ 



HL 

1.573 X 10“® 

1.474 X 10““ 

10.68 

1.601 X 10““ 

1.69 X 10““ 

9.47 

HV 

7.713 X 10“® 

1.139 X 10“® 

0.68 

-6.738 X 10““ 

1.24 X 10“® 

-0.54 

LV 

5.854 X 10“® 

1.042 X 10“® 

0.56 

1.131 X 10““ 

1.24 X 10“® 

0.09 

Combined 

1.541 X 10““ 

1.447 X 10““ 

10.65 

1.533 X 10““ 

1.659 X 10““ 

9.24 


Set 3 (BBH): 

= 6.975 X 10““ 



HL 

5.175 X 10“® 

1.474 X 10“® 

3.51 

4.725 X 10““ 

1.683 X 10““ 

2.81 

HV 

4.257 X 10“® 

1.139 X 10“® 

0.37 

-6.117 X 10““ 

1.407 X 10““ 

-0.43 

LV 

1.763 X 10“® 

1.042 X 10“® 

0.17 

3.968 X 10““ 

1.379 X 10““ 

-0.29 

Combined 

5.094 X 10““ 

1.447 X 10““ 

3.52 

4.448 X 10““ 

1.659 X 10““ 

2.68 


Set 4: (BNS-fBBH): 

Qc = 2.022 X 10““ 



HL 

2.056 X 10“® 

1.474 X 10““ 

13.94 

1.991 X 10““ 

1.684 X 10““ 

11.83 

HV 

9.674 X 10“'^ 

1.139 X 10“® 

0.85 

6.452 X 10““ 

1.352 X 10““ 

0.48 

LV 

7.211 X 10“® 

1.042 X 10“® 

0.69 

1.777 X 10““ 

1.328 X 10““ 

1.34 

Combined 

2.012 X 10““ 

1.448 X 10““ 

13.9 

1.968 X 10““ 

1.658 X 10““ 

11.87 


Observing scenario: 

0.0, = 2.756 X 10““ 



HL 

3.581 X 10“® 

6.869 X 10“"“ 

5.21 



— 

HV 

9.413 X 10“"“ 

6.207 X 10“® 

0.15 



— 

LV 

4.235 X 10“"“ 

5.723 X 10“® 

0.07 



— 

Combined 

3.505 X 10““ 

6.779 X 10“"“ 

5.17 

- - - 


test how much bias is added to the measurements of fla 
when we neglect to remove the loud detectable signals 
from the data streams^ compared to when we only use 
sub-threshold signals. By keeping the detectable events, 
we also increase as indicated in Table [IIT| We find 


^ In this MDSC we can choose to simply not include any indi¬ 
vidually detectable signals within the data streams. In reality, 
removing a detected signals from the data streams is very difficult 
as there may be some inaccuracy in measuring its true parameter 
used to produce the waveform which would leave some residual 
signal. Instead, we simply do not analyse the data that is known 
to contain the signal in the frequency band that we are searching. 


that in both cases we are able to measure the background 
estimate to within Icr of the true value, as well as obtain¬ 
ing similar SNRs, both greater than 10. 

With data set 3, where we consider a population of 
only BBH with a lower merger rate than before, we find 
that we are still able to make a detection but with a lower 
SNR of 3.52. This also gives the largest error in the mea¬ 
surement of fla with a measured value 1.3cr away from 
the true value. These results show that the estimation 
of fta is still possible given a non-continuous GW signal 
(see middle plot of Fig. [^. 

Data set 4, which is the combination of data set 1 and 
data set 3, gives a measured spectrum equal to the 
sum of the results from data sets 1 and 3. A plot of the 
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results for the three detector pairs, as well as the com¬ 
bined results, are shown in the left-hand plot of Fig. 
where the measured values flcare plotted in blue along 
with their respective error bars. Here the red dashed line 
shows the true value of Hq . We see here that we are able 
to recover the value of Hq to within 1.6% with the HL 
detector pair and 0.5% when we consider the combined 
results. This is likely due to the very high SNR of 13.9. 


B. Re-coloured 


The re-coloured data consists of the re-coloured initial 
LIGO and Virgo detector noise combined with the same 
data sets signal as for the Gaussian case. This anal¬ 
ysis more closely simulates the likely output from the 
advanced detectors, which will suffer from various envi¬ 
ronmental noise sources contaminating the data. The re¬ 
sults from the analysis of all the re-coloured data sets are 
presented in the right-hand section of Table III (next to 


the Gaussian results for easy comparison). As in the case 
of the Gaussian data set, set 0 contains only noise and 
serves as a baseline for the analysis. It has a combined 
point estimate well within Icr of 0. The results for data 
sets 1 and 2 are also consistent with the Gaussian sets. 
The SNR for these sets are about 20% lower than that of 
the Gaussian set, which is due to the non-Gaussian noise 
transients. They are also within Icr of the true values. 
Data set 3 and 4 show similar effects to data set 2, with 
a lower SNR than in the Gaussian case and with a recov¬ 
ered value of ila within 1.6 a of the true values for data 
set 3 and within Icr for data set 4. 


C. Observing Scenario 


We report the results from the observing scenario 
where we consider evolution of the detector sensitivities 
during different observational runs, the results of which 
are shown in Fig. In the left-hand plot we show the 
measured value of Hq, as a function of observation time 
with error bars included. The black dashed lines rep¬ 
resent the start/end of each phase as detailed in Sec¬ 
tion m and the red dashed line represents the true 
injected Dq, value while the solid black horizontal line 
represents Dq, = 0. In the right hand plot we show the 
combined SNR as a function of time for the same points 
in the left-hand plot, as well as the theoretical SNR given 
by Eq. (22) using the blue dashed line. Here again the 
vertical black dashed lines show the start/end of each 
phase and the red dashed line indicates SNR = 3 which 
we use as the threshold value for claiming a detection of 


the SGWB. 

The first resnlt to note is that over the course of the 
whole 5.5 years of the observing run the size of the er¬ 
ror bars reduces significantly. By the end of the first 1.5 
years of observations they have already reduced by over 
an order of magnitude. The second result to note is that 


already after the first 1.5 years of observation, we will ob¬ 
serve disagreement with the null result at 95% confidence 
(2a). The third result, which is in agreement with the 
theoretical model, is that we may be able to confirm the 
detection of the SGWB with an SNR > 3 after a period 
of about 3.5 years. At the end of the 5.5 year observing 
run, for this coalescence rate, we report that we have a 
total SNR of 5.15. We note that, in the right-hand plot 
of Fig. 1^ the measured SNR (0.95) for the end of the 
second phase is well above what is predicted by the theo¬ 
retical model (SNR = 0.35), although it is still within the 
Icr range. This is explained by the larger than average 
measurement of Hq, at the end of the second phase, as 
shown in the left-hand plot of Fig. which, whilst be¬ 
ing significantly larger than other measurements, is still 
within Icr of the true value. 


D. Parameter Estimation 


In order to construct parameter posterior distribu¬ 
tions, we use models for ^M(fi',9) for various sets of 
waveform parameters. We use a power-law and compact 
binary model for Eq- (24) is evaluated repeat¬ 

edly for each set of parameters and is maximized for those 
parameters that best fit the data. Parameter posterior 
distributions are constructed for parameter sets of equal 
likelihood. Example posteriors are shown in Fig. for 
the power-law and GBC models. Table [TV] shows results 
for all of the injection sets. We provide parameter esti¬ 
mates for the power-law model where, for CBG systems, 
the power-law index is a = 2/3. We also provide param¬ 
eter estimates and the true values for the CBG model. 
The constraints on the GBC background are relatively 
weak and highly dependent on the mass of the system 
(see the right hand plot of Fig. [^. Therefore, the limits 
we place on chirp mass and coalescence rate are in terms 
of bounds on the parameters. The bounds we place are 
consistent with the injected parameter values. We tested 
the case where we consider multiple CBG models, as was 
used to produce data set 4, and found that the posteriors 
are broadened by a significant amount. 


VI. CONCLUSION 


In this SGWB MDSC for the advanced detectors, 
LIGO and Virgo, we have presented our methods for 
the productions and analysis of multiple mock data sets, 
as well as the results and their scientific interpretation. 
We consistently find that the best results are obtained 
with the aLIGO detector pair (HL), compared to the 
aLIGO-AdVirgo pairs (HV, LV), though we still see some 
slight improvement when we consider the combined re¬ 
sults from all detector pairs. This is in agreement with 
what we expect given the difference in sensitivities and 
orientations of the pairs. In the case of onr three detector 
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FIG. 5: Left — Results for the year-long data set 4 with Gaussian noise where we show the measured value of with error 
bars included for each of the detector pairs as well as the combined result of the weighted sum of the three detector pairs. 
The horizontal red dashed line shows the true injected value of Qa = 2.022 x 10“® for this data set. Right — Results for the 
year-long set 4 with re-coloured noise where we show the measured value of Qa with error bars included for each of the detector 
pairs as well as the combined result of the weighted sum of the three detector pairs. The horizontal red dashed line shows the 
true injected value of fla = 2.022 x 10“® for this data set. 




Time (years) 


FIG. 6: Left — Gombined results for the estimated value of Lla as a function of time for the 5.5 long year observing scenario. The 
black dashed lines show the end point of each observing run, the red dashed line is the true injected value of Hq, = 2.756 x 10“®, 
the horizontal black solid line is Lie — 0 and the blue points are the measured Qa values with their error bars. Right — 
Gombined results for the measured SNR as a function of time for the 5.5 year long observing scenario. The black dashed lines 
show the end point of each observing run, the red dashed line represents SNR = 3, which we use as a threshold for claiming a 
detection, the blue dashed line shows the theoretical SNR as a function of time given by Eq. (221 and the blue points are the 
combined SNRs corresponding la error bars. 


network, the two aLIGO detectors are the best aligned 
and have the smallest separation. 

We have made comparisons between the use of Gaus¬ 
sian detector noise, which can be considered an idealistic 
case, and re-coloured noise data which is expected to be 
more realistic. We injected the same sources into both 
sets of noise to ensure that we are measuring the same 
signal. In both cases, we have been able to recover the 
injected value Hq to within 2a using each data set. 

From the analysis of the observing scenario data set 
we have shown that for the optimistic values of the GBG 


event rate the hrst deviation from the null result (at 95% 
confidence) may be observed as early as 1.5 years into 
the observation time. This assumes that we are able to 
achieve the designed sensitivities at the end of each ob¬ 
serving phase, as outlined in Section |IIB[ and that the 
coalescence rate of GBG signals is significantly large as 
to make it detectable within a few years of operations at 
design sensitivity. 

The results from data sets 1 and 3 have also shown 
that the theoretical prediction given in [55] holds true 
when applied to mock data. That is, the the statistical 
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TABLE IV: Parameter estimation results for the various data sets. We provide the 99% confidence limits for both the power- 
law and CBC models, as well as the injected parameters. The first column indicates the data set. The second column is the 
estimated amplitude of fla. The third column is the estimated power-law of the signal. The fourth and fifth columns gives the 
injected values of the amplitude and power-law. The sixth column give the estimated average chirp mass. The seventh column 
gives the estimated rate of events. The eighth and ninth columns gives the injected values of the average chirp mass and rate 
of events. 


Data Set 

^Ct 

a 

True Lie 

True a 

M 

A 

True M 

True A 

Set 1 

[4.0 X 10“^“, 4.0 X 10"**] 

[-0.3, 1.6] 

1.5 X 10“® 

0.65 

< 23 

> 1.0 

1.22 

10 

Set 2 

[4.0 X 10~®, 4.0 X 10“®] 

[-0.3, 1.6] 

1.5 X 10“® 

0.65 

< 22 

> 0.9 

1.22 

10 

Set 3 

[1.0 X 10“l^ 3.4 X 10“®] 

H, 2] 

7.6 X 10“^ 

0.73 

< 100 

> 0.04 

8.7 

0.3 

Set 4 

[8.7 X 10“®, 4.8 X 10“®] 

[0, 1.4] 

2.1 X 10“® 

0.68 

< 19 

> 1.6 

1.37 

10.2 




FIG. 7: Left — The posteriors for amplitude, Qa, and its spectral index, a, for data set 1 with the 99% confidence level (blue), 
95% confidence level (red) and 68% confidence level (black) shown. We show the true parameter values with an “X”. Right — 
The posteriors for average chirp mass, A4, and coalescence rate. A, for data set 1 with the 99% confidence level (blue), 95% 
confidence level (red) and 68% confidence level (black) shown. We denote the correct CBC parameters by an “X”. 


properties of the CBC GW signals, whether it be a con¬ 
tinuous signal or more popcorn like, do not matter when 
we make a measurement of all that is important is 
the total number of events that coalesce within the ob¬ 
servational period and the GW energy spectrum emitted 
by each event. 

Finally, we have shown that we are able to apply pa¬ 
rameter estimation methods to the data in order to place 
confidence levels on different parameters. The detection 
of a stochastic signal will not be able to provide enough 
information by itself to place tight constraints on these 
parameters, but when considered in combination with 
detections of single events, it can become a very useful 
tool to explore the ensemble of sources from the whole 
universe |44) . 

The results from the estimation of the average chirp 
mass and coalescence rate also show that there is equal 
probability of having a high rate of events and low av¬ 
erage mass as having a low rate of events and high av¬ 
erage mass. In both cases the amplitude of the signal 
and and the spectral index will be the same for the fre¬ 
quency range we search over, but in the first case the 


GW signals will be continuous while in the second case 
the GW signals will be highly non-Gaussian and popcorn 
like. The isotropic CC search we implement here is in¬ 
sensitive to two types of signals as it just considers the 
average strength of the signal over the observational pe¬ 
riod. In order to be able to differentiate these two signal 
types a non-Gaussian analysis must be developed that is 
able to search over both time and frequency [73] . 

Future MDSCs may wish to explore several areas that 
have not been covered here, such as, the inclusion of inter¬ 
mediate mass black holes (IMBH), which may coalesced 
in the middle of the frequency search band, and there¬ 
fore, given a high enough rate, may affect the analysis. 
Or, we could add a loud SGWB signal of cosmological ori¬ 
gin that has a spectral index that differs from that of the 
astrophysical contribution. Another important question 
would be to investigate the behaviour of the CC analysis 
when correlated noise between different pairs of detectors 
is included in the mock data |65L I66j . The continuation 
of MDSCs will be an important part of the verification 
process for LIGO and Virgo when a SGWB is eventually 
observed. 





















14 


Acknowledgments 

DM acknowledges the PhD financial support from the 
Observatoire de la Cote dAzur and the PACA region. 
MC was supported by the National Science Founda¬ 
tion Graduate Research Fellowship Program, under NSF 
grant number DGE 1144152. SM and JDR were sup¬ 
ported NSF grants GREST HRD-1242090 and PHY- 
1205585. JDR also acknowledges support from NSF 


grants PHY-1505861 and PHY-1066293, and the hos¬ 
pitality of the Aspen Center for Physics. NEC is sup¬ 
ported by NSF grants PHY-1505373 and PHY-1204371. 
VM is supported in part by NSF grant PHY-1204944 
at the University of Minnesota. The stochastic GW 
search has been carried out using the MatApps soft¬ 
ware available at https://www.lsc-group.phys.uwm. 
edu/daswg/proj e ct s/mat apps.html 


[1] L. P. Grishchuk, Sov. Phys. JETP 40, 409 (1975). 

[2] A. A. Starobinskii, JETP Lett. 30 , 682 (1979). 

[3] L. P. Grishchuk, Phys. Rev. D 48, 3513 (1993). 

[4] M. Gasperini and G. Veneziano, Astropart. Phys. 1 , 317 
(1993). 

[5] A. Buonanno, M. Maggiore, and G. Ungarelli, Phys. Rev. 
D 55, 3330 (1997). 

[6] J. F. Dufaux, D. G. Figueroa, and J. Garcia-Bellido, 
Phys. Rev. D 82, 083518 (2010). 

[7] T. Damour and A. Vilenkin, Phys. Rev. D 71 , 063510 
(2005). 

[8] X. Siemens, V. Mandic, and J. Creighton, Phys. Rev. 
Lett. 98, 111101 (2007). 

[9] S. Olmez, V. Mandic, and X. Siemens, Phys. Rev. D 81 , 
104028 (2010). 

[10] T. Regimbau, S. Giampanis, X. Siemens, and V. Mandic, 
Phys. Rev. D 85, 066001 (2012). 

[11] C. Caprini, R. Durrer, and G. Servant, Phys. Rev. D 77 , 
124015 (2008). 

[12] C. Caprini, R. Durrer, T. Konstandin, and G. Servant, 
Phys. Rev. D 79 , 083519 (2009). 

[13] C. Caprini, R. Durrer, and G. Servant, Journal of Cos¬ 
mology and Astroparticle Physics 12 , 024 (2009). 

[14] M. Maggiore, Physics Reports 331 , 283 (2000). 

[15] P. Binetruy, A. Bohe, C. Caprini, and J.-F. Dufaux, Jour¬ 
nal of Cosmology and Astroparticle Physics 6, 27 (2012). 

[16] T. Regimbau, Res. Astron. Astrophys. 11 , 369 (2011). 

[17] A. Buonanno, G. Sigl, G. G. Raffelt, H. T. Janka, and 
E. Muller, Phys. Rev. D 72 , 084001 (2005). 

[18] P. Sandick, K. A. Olive, F. Daigne, and E. Vangioni, 
Phys. Rev. D 73 , 104024 (2006). 

[19] S. Marassi, R. Schneider, and V. Ferrari, MNRAS 398 , 
293 (2009). 

[20] X. J. Zhu, E. Howell, and D. Blair, MNRASL 409, L132 

( 2010 ). 

[21] T. Regimbau and J. A. de Freitas Pacheco, Astron. & 
Astrophys. 376 , 381 (2001). 

[22] P. A. Rosado, Phys. Rev. D 86, 104007 (2012). 

[23] T. Regimbau and J. A. de Freitas Pacheco, Astron. & 
Astrophys. 447 , 1 (2006). 

[24] E. Howell, T. Regimbau, A. Corsi, D. Coward, and 
R. Burman, MNRAS 410 , 2123 (2011). 

[25] S. Marassi, R. Ciolfi, R. Schneider, L. Stella, and V. Fer¬ 
rari, MNRAS 411 , 2549 (2011). 

[26] C. J. Wu, V. Mandic, and T. Regimbau, Phys. Rev. D 
87 , 042002 (2013). 

[27] J. C. N. de Araujo and G. F. Marranghello, General Rel¬ 
ativity and Gravitation 41 , 1389 (2009). 

[28] V. Ferrari, S. Matarrese, and R. Schneider, MNRAS 303, 


258 (1999). 

[29] E. Howell, D. Coward, R. Burman, D. Blair, and 

J. Gilmore, MNRAS 351, 1237 (2004). 

[30] X. J. Zhu, X. L. Fan, and Z. H. Zhu, Astrophys. J. 729 , 
59 (2011). 

[31] X. J. Zhu, E. Howell, T. Regimbau, D. Blair, and Z. H. 
Zhu, Astrophys. J. 739 , 86 (2011). 

[32] P. A. Rosado, Phys. Rev. D 84, 084004 (2011). 

[33] S. Marassi, R. Schneider, G. Corvino, V. Ferrari, and 
S. P. Zwart, Phys. Rev. D 84, 124037 (2011). 

[34] C. Wu, V. Mandic, and T. Regimbau, Phys. Rev. D 85, 
104024 (2012). 

[35] X.-J. Zhu, E. J. Howell, D. G. Blair, and Z.-H. Zhu, MN¬ 
RAS 431 , 882 (2013). 

[36] J. Aasi et al. (LIGO Scientific Collaboration), Class. 
Quantum Grav. 32 , 074001 (2015). 

[37] F. Acernese et al. (Virgo Collaboration), Class. Quantum 
Grav. 32 , 024001 (2015). 

[38] B. P. Abbott et al. (LIGO Scientific Collaboration), Rep. 
Prog. Phys. 72 , 076901 (2009). 

[39] F. Acernese et al. (Virgo Collaboration), AIP Conf. Proc. 
794 , 307 (2005). 

[40] B. P. Abbott et al. (LIGO Scientific Collaboration and 
Virgo Collaboration), Nature 460 , 990 (2009). 

[41] T. Regimbau and S. A. Hughes, Phys. Rev. D 79 , 062002 
(2009). 

[42] E. Thrane, S. Ballmer, J. D. Romano, S. Mitra, 
D. Talukder, S. Bose, and V. Mandic, Phys. Rev. D 80 , 
122002 (2009). 

[43] I. Kowalska, T. Regimbau, T. Bulik, M. Dominik, and 

K. Belczynski, Astron. & Astrophys. 574, A58 (2015). 

[44] V. Mandic, E. Thrane, S. Giampanis, and T. Regimbau, 
Phys. Rev. Lett. 109 , 171102 (2012). 

[45] J. Abadie et al. (LIGO Scientific Collaboration and 
Virgo Collaboration), Class. Quantum Grav. 27 , 173001 
( 2010 ). 

[46] T. Regimbau, T. Dent, W. Del Pozzo, S. Giampa¬ 
nis, T. G. F. Li, C. Robinson, C. Van Den Broeck, 
D. Meacher, C. Rodriguez, B. S. Sathyaprakash, et al., 
Phys. Rev. D 86, 122001 (2012). 

]47] T. Regimbau, D. Meacher, and M. Coughlin, Phys. Rev. 
D 89, 084046 (2014). 

[48] J. Aasi et al. (LIGO Scientific Collaboration and Virgo 
Collaboration), Phys. Rev. Lett. 113, 231101 (2014). 

[49] J. Abadie et al. (LIGO Scientific Collaboration and Virgo 
Collaboration), Phys. Rev. D 82 , 102001 (2010). 

[50] B. P. Abbott et al. (LIGO Scientific Collaboration), 
Phys. Rev. D 79, 122001 (2009). 

[51] J. Aasi et al. (LIGO Scientific Collaboration and Virgo 



15 


Collaboration), arXiv 1304, 0670 (2013). 

[52] T. Regimbau, K. Siellez, D. Meacher, B. Gendre, and 

M. Boer, Astrophys. J. 799, 69 (2015). 

[53] K. Belczynski, R. Perna, T. Bulik, V. Kalogera, 

N. Ivanova, and D. Q. Lamb, Astrophys. J. 648, 1110 
(2006). 

[54] E. Berger and other, Astrophys. J. 664, 1000 (2007). 

[55] K. Belczynski and V. Kalogera, Astrophys. J. Lett. 550, 
L183 (2001). 

[56] T. Bulik, K. Belczynski, and B. Rudak, Astron. & As¬ 
trophys. 415, 407 (2004). 

[57] A. M. Hopkins and J. F. Beacom, Astrophys. J. 651, 142 
(2006). 

[58] K. Belczynski, V. Kalogera, F. A. Rasio, R. E. Taam, 
A. Zezas, T. Bulik, T. J. Maccarone, and N. Ivanova, 
Astrophys. J. Suppl. 174, 223 (2008). 

[59] M. Dominik, K. Belczynski, C. Fryer, D. E. Holz, 
E. Berti, T. Bulik, 1. Mandel, and R. O’Shaughnessy, 
Astrophys. J. 759, 52 (2012). 

[60] S. Oslowski, T. Bulik, D. Gondek-Rosinska, and K. Bel¬ 
czynski, MNRAS 413, 461 (2011). 

[61] B. P. Abbott et al. (LIGO Scientific Collaboration), 
Phys. Rev. D 80, 047101 (2009). 


[62] J. Abadie et al. (LIGO Scientific Gollaboration and Virgo 
Gollaboration), Phys. Rev. D 85, 082002 (2012). 

[63] P. Ajith et al., Phys. Rev. D 77, 104017 (2008). 

[64] A. Buonanno, B. R. Iyer, E. Ochsner, Y. Pan, and B. S. 
Sathyaprakash, Phys. Rev. D 80, 084043 (2009). 

[65] E. Thrane, N. Ghristensen, and R. M. S. Schofield, Phys. 
Rev. D 87, 123009 (2013). 

[66] E. Thrane, N. Christensen, R. M. S. Schofield, and A. Ef- 
fler, Phys. Rev. D 90, 023013 (2014). 

[67] D. Coward and T. Regimbau, New Astron. Rev. 50, 461 
(2006). 

[68] D. Meacher, E. Thrane, and T. Regimbau, Phys. Rev. D 
89, 084063 (2014). 

[69] B. Allen and J. D. Romano, Phys. Rev. D 59, 102001 
(1999). 

[70] N. Christensen, Phys. Rev. D 46, 5250 (1992). 

[71] B. Allen, in Relativistic Gravitation and Gravitational 
Radiation, edited by J. A. Marck and J. P. Lasota (1997), 
p. 373. 

[72] J. Abadie et al. (LIGO Scientific Gollaboration and Virgo 
Collaboration), Phys. Rev. D 85, 122001 (2012). 

[73] E. Thrane, Phys. Rev. D 87, 043009 (2013). 



